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Abstract 

The Gowdy spacetimes are vacuum solutions of Einstein's equations with two 
commuting Killing vectors having compact spacelike orbits with T 3 , S 2 x S 1 or S 3 
topology. In the case of T 3 topology, Kichenassamy and Rendall have found a family 
of singular solutions which are asymptotically velocity dominated by construction. 
In the case when the velocity is between zero and one, the solutions depend on the 
maximal number of free functions. We consider the similar case with S 2 x S 1 or 
S 3 topology, where the main complication is the presence of symmetry axes. We 
use Fuchsian techniques to show the existence of singular solutions similar to the 
T 3 case. We first solve the analytic case and then generalise to the smooth case 
by approximating smooth data with a sequence of analytic data. However, for the 
metric to be smooth at the axes, the velocity must be 1 or 3 there, which is outside the 
range where the constructed solutions depend on the full number of free functions. 
A plausible explanation is that in general a spiky feature may develop at the axis, a 
situation which is unsuitable for a direct treatment by Fuchsian methods. 



1 Introduction 

The singularity theorems by Hawking and Penrose show that in general spacetimes, singu- 
larities will be present [31]. Apart from existence, these results give very little information 



on the nature of these singularities. 

There is however a general proposal for the structure of the generic cosmological sin- 
gularity, put forth by Belinskii, Khalatnikov and Lifshitz (BKL) The argument is 
based on formal calculations which originally were without rigorous mathematical proofs. 
In essence, the BKL picture consists of two predictions: firstly, that the evolution of a 
spatial point near the singularity is unaffected by the evolution of nearby points. Thus the 
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dynamics of an inhomogeneous model can be described by a homogeneous model at each 
spatial point, with the parameters of the homogeneous model depending on the spatial 
coordinates. Secondly, the dynamics are expected be oscillatory, similar to the situation 
in the homogeneous Bianchi IX models. Loosely speaking, the evolution of a spatial point 
near a cosmological singularity should appear as an infinite sequence of epochs, each similar 
to a different Kasner spacetime. 

Some cosmological models exhibit a simpler behaviour without oscillations, termed 
'asymptotically velocity dominated' (AVD) |TJ, 25|. For these models there is no oscillation 



and the evolution of a given spatial point approaches a particular Kasner solution. This is 
exhibited in the evolution equations by the spatial derivative terms becoming negligible in 
relation to time derivatives. 

The behaviour of general spacetimes is still out of reach of the mathematical tech- 
niques available, so it makes sense to study particular subclasses of solutions, subject to 
some particular symmetry and/or a particular choice of matter. Recently there has been 



considerable progress in the study of spatially homogeneous spacetimes |3q, B7l 40, 41, 43 



There is also an interesting result which establishes AVD behaviour for the general Einstein 
equations coupled to a scalar field, without any symmetry assumptions 0. 

It is of obvious interest to relax the symmetry assumptions from homogeneous to inho- 
mogeneous spacetimes. For simplicity, we restrict attention to the spatially compact case 
in vacuum. The largest symmetry group that is compatible with inhomogeneities is then 
{7(1) x U(l), which in the spatially compact case leads to the Gowdy models |L7], [HJ. The 
possible spatial topologies are then essentially T 3 , S 2 x S 1 or S 3 . (In the more general case 
of local {7(1) x £/(!) symmetry, other topologies are possible, see 



For Gowdy spacetimes it is possible to define a quantity, the 'velocity', which describes 
the dynamics near the singularity. The terminology is slightly misleading since the velocity 
is only uniquely defined up to a sign. While a sign can be introduced, it depends on 
the particular parametrisation of the metric. If AVD holds, the velocity has a (space 
dependent) limit at the singularity. 

For the subclass of polarised Gowdy models, which describe universa with polarised 
gravitational waves, the Einstein evolution equations reduce to a single linear ordinary 
differential equation, which can be solved formally in terms of a series of Legendre poly- 
nomials. The polarised solutions have been shown to be AVD fll2| , p5| , p0| . 



For the full Gowdy models with T 3 spatial topology, numerical simulations show that 
the velocity has a limit between and 1 for almost all spatial points in generic situations 
H H, |22|, p^] . The exceptions are isolated points where the velocity remains outside this 
range, which results in a discontinuous behaviour, 'spikes', in the limit. A family of singular 
solutions have been constructed analytically by Kichenassamy and Rendall using Fuchsian 



methods ||28|| , based on previous work on formal asymptotic expansions by Grubisic and 
Moncrief These solutions have a prescribed asymptotic behaviour. In the case when 
the asymptotic velocity is between and 1, the solutions depend on the full number of free 
functions, which indicates that they correspond to an open set of initial data. Recently, 
a similar family of solutions with spikes have been constructed by Rendall and Weaver, 



showing that the spikes are real features of the model and not numerical artefacts |5U 
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We should also mention that there is a result on the existence of Gowdy T 3 solutions with 
small data close to that of a Kasner (2/3,2/3, —1/3) spacetime 

For the unpolarised Gowdy spacetimes with S 2 x S 1 or S 3 topology, much less is known. 
Garfmkle has done numerical simulations of Gowdy S 2 x S 1 spacetimes, which show similar 
behaviour as in the T 3 case, including the appearance of spikes. It is the purpose of this 
paper to provide a Fuchsian analysis of Gowdy S 2 x S 1 and S 3 spacetimes similar to that 



in 28 



for T 3 spacetimes. 

The Fuchsian algorithm [27, 28 1 has proven to be a valuable tool in showing the existence 
of singular solutions to nonlinear wave equations. Variants of it has been used to show 
the existence of Cauchy horizons with certain properties | 30| , to analyse the properties of 
isotropic singularities |33|, |34|, [13], ||, f|, || , to construct families of singular Gowdy or plane 
symmetric spacetimes p8| , ^, |5], [32], |23| and for the general result on scalar field spacetimes 
I- 

The original version of the Fuchsian technique only applies to the analytic case, since 
the proof uses the Cauchy formula in the complex domain to estimate derivatives in terms 
of functions themselves. There are a number of results extending the Fuchsian algorithm 

[38] in various contexts. In section |] below we 



to give smooth solutions []33], |34|, [L3|, [26|. 
will modify the argument of [3S] to apply to our case. 

Finally, we note that there is an additional motivation for studying Gowdy spacetimes 
with S 2 x S 1 topology apart from cosmology, since axisymmetric and stationary black hole 
interiors are in fact Gowdy S 2 x S 1 spacetimes. These include the region of Kerr spacetime 
between the inner and outer horizon |35|] and the interiors of the 'distorted black holes' 

m. 



studied by Geroch and Hartle 

The outline of the paper is as follows. In section |2] we give an introduction to the 
Gowdy spacetimes, introduce coordinates and a parametrisation of the metric, and rewrite 
the field equations in a suitable form. We also discuss the restrictions needed for the 
metric to be smooth at the axes and the properties of some previously known solutions. In 
section |^ we give a short review of the Fuchsian technique of ||28|| , and construct analytic 
solutions of the field equations. This is then extended to smooth solutions in section f|, by 
means of a generalisation of the technique in [35|. Finally we discuss the shortcomings and 
implications of the results in section |5|. 



2 Gowdy spacetimes 

We assume that {M,g) is a U(l) x ?7(l)-symmetric spatially compact spacetime with 
spacelike group orbits. More precisely, let M = I x E where / C R and E is a connected 
orientable compact three manifold, and assume that the U(l) x U(l) group acts effectively 
on E. It then follows [TI| that the topology of E is that of either T 3 , S 2 x S 1 , S 3 or a lens 
space L(p, q). Since the T 3 case is treated in [^8] and ]38| and L(p,q) may be covered by 
S 3 , we only consider E ~ S 2 x S 1 or S 3 here.[] 

1 J28[ and [Q treat the T 3 case with vanishing twist. See |24| for a similar treatment with nonvanishing 
twist but restricted to the polarised case. 
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We introduce coordinates on E as follows. Let 9 label the orbits of U(l) x U(l) in 
£ and let (0, x) be coordinates on each orbit induced by the standard coordinates on 
U(l) x U(l). For S ~ S 2 x S 1 , we choose x to be a cyclic coordinate on the S 1 part, 
X G [0, 27r] mod 2tt, and (6 1 , 0) to be spherical coordinates on the S 2 part, 9 G [0, tt] and 
4> G [0, 27r] mod 27r. This fixes (0, x) up to translation. For £ ~ S 3 , the fix points to the 
first U(l) factor of U(l) x £7(1) form a circle Si, and the set of fix points to the second 
factor is a similar circle S2 (this is a nontrivial consequence of the effective action, see, 
e.g., PI]). The parametrisation is chosen such that Si and S 2 correspond to 9 = and 



9 = tt, respectively. This fixes ((f), x) up to translation in this well. Note that the 

two cases are similar in a neighbourhood of one of the axes. 

If (M, g) is assumed to be a maximal globally hyperbolic development of generic data 
on S, then M contains the set (0, tt) x E, on which the line element of the metric g may 
be parametrised as 

ds 2 = e A (-dt 2 + d9 2 ) + da 2 , (1) 

where da 2 is a metric on the orbits with determinant R = csintsin# for some positive 
constant c [[TI|. We set c = 1 from now on, since a constant conformal factor does not affect 
the qualitative aspects of the model. Here 'generic' refers to an open and dense subset of 
all Cauchy data; see JTO] for the exact definition. We will call this model the Gowdy model 



over 



0- 



In the T 3 case, the metric can also be written in the form ([!]), but with orbit metric 
determinant t. It follows that for S 2 x S 1 or S 3 , the metric can be brought into the T 3 



form by a change of coordinates whenever 9 7^ and 9 ^ tt. Thus the result of and 



| 38fl may be applied locally on any subset of M where 9 is bounded away from and tt, and 
what remains is to study a neighbourhood of one of the symmetry axes. (Strictly speaking, 
in the smooth case it is necessary to show that local solutions can be patched together. 
This may be done using the domain of dependence result of Theorem |4.1| below.) Also, 
by the symmetry of the problem, we may restrict attention to a neighbourhood of 9 = 
and t = 0. When we write 'the axis' below, we will mean the axis at 9 = 0, although all 
arguments apply to the other axis at 9 = tt as well. 

2.1 Parametrisations 

The orbit metric da 2 may be parametrised in several ways. If we use the parametrisation 

ds 2 = e A (-dt 2 + d9 2 ) + R[e p (d<P + Qd X ) 2 + e~ p d X 2 ] , (2) 

the Einstein equations decouple as two evolution equations involving P and Q and two 
constraint equations involving P, Q and A |10| . Let (t, 9) G M C R 2 with the Minkowski 



metric 77 = —dt 2 + d9 2 and (P, Q) E N C R 2 with the hyperbolic metric h = dP 2 + e 2P dQ 2 . 
Given the map f:M9(^)H (P, Q) e N, we define a section d^ of T*M ® ^*(TN) by 

d^ = dP®^*— + dQ®^*— . (3) 
oP oQ 
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Then the evolution equations may be written as the wave-map type equation 



tr„(D(M)(i2d*)) =0, (4) 

where T>(h,rj) is a connection on T*M ® ty*(TN) defined in terms of the connections of 
(M,rj) and (N, h) (see pf for details). 

The 'velocity' u(t, 6q) is defined up to a sign as the ^-velocity of the solution (P, Q) of 
(|) along a curve 9 = 9q, rescaled by an appropriate factor. In the T 3 case, this factor is 
t, but in our case it is more convenient to take the factor to be tant. So we define the 
velocity to be 



v 



^(DPy + e* p (DQ)\ (5) 



where D = tant d/dt. 

An 'asymptotically velocity dominated' (AVD) solution is a solution which asymp- 
totically approach a solution to the 'AVD equations', obtained by dropping pure spatial 



derivative terms in the evolution equations [14, E5[. For an AVD solution, it follows im- 



mediately from the AVD equations that the velocity v tends to a finite limit as t — > 0. We 
will refer to this limit as the asymptotic velocity. It is also possible to fix the sign of the 
velocity to be the same as the sign of DP, note however that this sign depends upon the 
parametrisation as we will see below. 

Note that P cannot be regular at the axis. The regular parametrisation used by Garfin- 
kle in his numerical work [15 is obtained from (0) by setting 



P = P + In sin 9 and A = P + 7 + In sin t, (6) 

with 7, P and Q smooth functions of t and 9. Unfortunately, the corresponding evolution 
equations contain coefficients that are singular at the axes. This is not a major obstacle 
for the numerics since the full terms appearing in the equations are well behaved, but it 
does prevent a direct analytical treatment. 

Another parametrisation may be obtained by interchanging the roles of the Killing 
vectors d/d<ft and d/dx- This corresponds to an inversion in the hyperbolic plane (N, h), 
so the field equations are invariant under this reparametrisation. Explicitly, 

ds 2 = e A (-dt 2 + d9 2 ) + R[e Y (d X + Xd<f)) 2 + e~ y d0 2 ] , (7) 

where 

Y = P + He~ 2P + Q 2 ) and X = e _J* + (8) 

Note that there is a subtlety here. If P — > 00 as t — > 0, the rate of blow-up of Y will be 
dramatically different depending on whether Q vanishes or not. This is a consequence of 
the parametrisation, and must be distinguished from real geometric features of the model. 
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The parametrisation (|7|) is again singular at the axis, but this may be resolved by 
substituting Y = Z — InR. We also set A = A — Z + 2 In sin t. The metric then becomes 

ds 2 = e x ~ z sin 2 1 (-dt 2 + d9 2 ) + e z (d X + Xd<f)) 2 + R 2 e~ z d<p 2 . (9) 

From Lemma 5.1 of |Kj, (§) defines a smooth metric with the desired symmetry if and 
only if A, Z and X are smooth functions of 9 on S 2 , with X and A vanishing at 9 = and 
9 = 7T. Note that a function is smooth in a neighbourhood of # = or 9 = tt in S 2 if and 
only if it is a smooth function of sin 2 9 there. 

Because of the invariance of h under inversion, the velocity is 



v = s/{DY) 2 + e 2Y (DX) 2 = y/(l - DZ) 2 + e 2Z R~ 2 {DX) 2 . (10) 

As mentioned above, it is possible to assign a sign to v by setting it equal to the sign of DP. 
Assume, for illustrative purposes, that Q is independent of t and that DP — > oo as t — > 0. 
Then it follows from (D that DZ -> DP as t -> wherever Q ^ 0, but £>Z -> -DP 
for those values of # where Q vanishes. So the sign defined with respect to DP cannot be 
expressed in terms of DZ alone. Also, if DP tends to a smooth function of x as t — > 0, the 
corresponding limit of DZ has discontinuities at points where Q has isolated zeros (note 
that for any smooth metric of the form (0), Q must vanish at the axes). Since DP has a 
perfectly regular asymptotic behaviour in this case, these discontinuities are artefacts of 



the parametrisation (called 'false spikes' in [39]). The same argument may of course also 



be applied in the other direction, giving discontinuities in the asymptotic behaviour of DP 
from a well-behaved DZ. For the rest of this paper, we fix the sign of the velocity to be 
equal to the sign of 1 — DZ. 

2.2 Einstein's equations 

In the last section we obtained a parametrisation which is regular at the axis. We will now 
modify the parametrisation to obtain evolution equations with regular coefficients as well. 

We will denote partial derivatives by subscripts, e.g., fg := dgf = df /d6. As mentioned 
above, Einstein's equations decompose into the wave map equation (H) for Z and X and 
two constraint equations for A, Z and X. If we denote the covariant derivative of the 
Minkowski metric rj on iV by T>, @ can be written as 

FC l V A (RV A Z) = R- 2 e 2Z V A XV A X, (11a) 
RV A (R- 1 e 2Z V A X) = 0, (lib) 

and the constraint equations are 

2(A± ± cot t)R± = 4(Ro)± + R(Zl + e 2Z R- 2 X 2 ± ), (12) 



where we have used the notation f± := f e ± f t [ID]. 

The integrability of the constraints ([12]) are ensured by the evolution equations ([TT|). 
When R + and R- are nonzero, the constraints determine A up to a constant, and the 
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constant is fixed by the smoothness requirement, i.e., that A = at 9 = 0. Since R = 
sin t sin 9, R + and R_ vanish when 9 = t and 9 = n — t respectively. So this happens for 
two values of 9 in each Cauchy surface. At these points, the constraints become 'matching 
conditions' for Z and X JT^, f42"fl . If the matching conditions hold on an initial Cauchy 
surface, they are preserved by the evolution equations In what follows we will assume 
without further comment that the solutions are chosen such that the matching conditions 
hold. 

To obtain a system with smooth coefficients, we reparametrise the metric using an Ernst 
potential as in [T(| , a procedure which is similar to the Kramer-Neugebauer transformation 
in stationary and axisymmetric spacetimes p9j . Let 

tt ■= -R- 1 e 2Z (X e dt + X t d9). (13) 



It then follows from (|1 lb|) that is a closed form. Thus we may write = duo for some 
function oo, defined up to a constant, on any simply connected open set containing 9 = 0. 
Inverting the relation (|T3| ) shows that if uo is a smooth function of 9 on a neighbourhood of 
9 = on S 2 , so is X. This determines X up to a constant, which is fixed by the requirement 
that X = at 9 = so that the metric is smooth at the axis. Following |39fl , we will call 
the transformation (Y,X) i— > (Z,uj) a 'Gowdy-to- Ernst' transformation. 
Expressing the evolution equation ( |lla| ) in terms of Z and u gives 

R~ 1 V A (RV A Z) = -e- 2Z V A uV A uj, (14a) 

and the identity T>\ A T) B -\X = together with (|i~3"D gives a second evolution equation 

V A (Re' 2Z V A uo) = 0. (14b) 

Expanding (|TJ) we get 

Z a + cot t Z t - Zee - cot 9 Z e = -e~ 2Z ' (uj 2 - coq), (15a) 
ujtt + cottu t — ujee — cot^cug = 1{Z t uj t — ZqOJq). (15b) 

We will now rewrite the equations in a form which makes the regularity of the coef- 
ficients at the axis explicit, and which is more suited for the Fuchsian techniques. Put 
t := sint and D := tantd t = rd T as before, and let A := d 2 . + cot9de + (sin^) -2 ^, 
e.g., the Laplacian on S 2 with respect to the metric d# 2 + sin 2 #d0 2 induced by the nat- 
ural embedding in Euclidean R 3 . We also write V/ ■ Vg := fege + (sin^) -2 and 
(V/) 2 := V/ • V/. Since Z and to are independent of 4>, (f[|) can be written as 

(1 - t 2 )D 2 Z - r 2 DZ = r 2 AZ - e~ 2Z [(1 - r 2 ) {Duo) 2 - t 2 {Vuj) 2 ] , (16a) 
(1 - t 2 )D 2 uj - t 2 Duj = t 2 Auj + 2 [(1 - t 2 )DZDuj -t 2 VZ- Vuo] . (16b) 

The constraints (O) expressed in terms of Z and u are 



DX + tan 2 1 cot 9 Xg = -2(1 + tan 2 1) + ^ (DZ) 2 + tan 2 1 Z 2 e + e~ 2Z {(Duo) 2 + tan 2 1 ufj 

(17a) 

DX + tan 9X e = tan 9 (Z e DZ + e" 2 ^ Duo) . (17b) 
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In the analytic case, the Fuchsian techniques can in fact be applied directly to fllED in 
the variables r and 9. For the smooth case however, we have to cast the equations into 
symmetric hyperbolic form, which cannot be done using the 9 variable alone because the 
axial symmetry will lead to coefficients singular in 9 at 9 = 0. We therefore introduce 
an approximately Cartesian coordinate system on a neighbourhood of 9 = on S 2 . Let 
x := sin 9 cos <fi and y := sin 9 sin <p. In the variables (x,y), the metric on S 2 is 

d9 2 + sin 2 9 d<p 2 = (1 - x 2 - y 2 )' 1 [(1 - y 2 )dx 2 + 2xy dxdy + (1 - x 2 )dy 2 ] , (18) 

hence 

A = (1 - x 2 ) d 2 x - 2xy d x d y + (1 - y 2 ) d 2 - 2x d x - 2y d y , (19) 

and 

V/ • V<7 = (1 - x 2 )f x g x - xy(f x g y + f y g x ) + (1 - y 2 )f y g y . (20) 

Inserting these expressions in (|i6|) then gives a system in the coordinates t, x and y with 
coefficients regular in x and y at the axis. 

2.3 Restrictions at the axes 

As mentioned in the introduction, asymptotic velocity dominance for Gowdy spacetimes 
may be loosely formulated as that the metric tends to a Kasner metric when t — > for 
a fixed point in E. The presence of symmetry axes imposes additional restrictions on the 
asymptotic behaviour along the axis. 

To make this precise we calculate the generalised Kasner exponents. These are defined 
as the eigenvalues g« of the renormalised second fundamental form (tr k)~ l kij, expressed in 
an orthonormal frame on E. From the expression for the metric (§) we find that 

q 1 = (DX- DZ + 2)/(DX- DZ + 4), (21a) 
q 2>3 = {l±v)/{DX-DZ + 4), (21b) 

where v is the velocity as given in ([10|) . 

Because of the rotational symmetry, the only possible values of the Kasner exponents 
at the axis are (0,0,1) and (2/3,2/3,-1/3) since the two eigenvalues corresponding to 
eigenvectors tangent to S must agree. From the smoothness requirements DX and DX 
vanish at the axis, and using our sign convention we have v = 1 — DZ there. It then follows 
immediately from (^l|) that for the (0, 0, 1) case, DZ — > 2 and v — > — 1 as t — * along the 
axis, while for the (2/3,2/3, —1/3) case, DZ — > —2 and v — > 3. 

As noted at the beginning of this section, the T 3 results of [^8j applies to sets not 
containing the axis at 9 = 0. In that case, we obtain solutions with the full number of 
free functions only in the case when < v < 1 (and DP < 0) in the limit t — > 0. Also, 
numerical simulations have shown that the velocity is indeed driven into this range, except 
at isolated spatial points (see also |39| for a heuristic argument). So we cannot hope 
to show the existence of velocity dominated solutions depending on the full number of free 
functions since v is either —1 or 3 on the axis. (Actually, in our case DZ will play the role 
of v, but the same conclusion can be drawn since v = 1 — DZ on the axis.) 



S 



2.4 Some known solutions 



2.4.1 Polarised solutions 

The polarised solutions form the subclass with Q = 0, so called because they describe 
spacetimes with polarised gravitational waves. In this case the evolution equations flToD 
reduce to a linear ordinary differential equation for P, the Euler-Poisson-Darboux equation. 
The solutions may be expressed explicitly in terms of Legendre polynomials [20|. Also, 



rigorous asymptotic expansions in a neighbourhood of the singularity have been found 



12 . In our notation, 



Z(t,6) = k{6)kismt + (p(6) +u{t,6), (22a) 
A(t,0) = ^{k 2 -4)lnsint + V(0) +v{t,6), (22b) 

with \D m dgu\ < C sin 2 i|lnsini| and \D m d^v\ < C sin 2 1 1 In sin 1 1 2 for m = 0, 1, all n and 
a constant C. The functions k and (p may be chosen freely as long as they satisfy the 
constraints, and the solutions are all asymptotically velocity dominated. Since A must 
vanish at the axis for the metric to be smooth, it follows immediately from the expansions 
that k = ±2 there, and the velocity u = 1 — k must tend to —1 or 3 along the axis. The 
Kretschmann curvature scalar is unbounded as t — > unless v = 1 — k — ±1, kg — and 
k ee = 0. 

There is a way in which one may obtained unpolarised solutions from a given polarised 
solution by means of an Ehlers transformation. This is used to obtain the 'reference 
solution' used to validate the numerical code in |T3| (see also [IB]). The technique works 



because a combination of an Ehlers transformation with Gowdy-to-Ernst transformations 
gives an isometry of the hyperbolic plane (N, h). The velocity of the transformed solution 
has the same asymptotic behaviour as for the original polarised solution. 

2.4.2 Black hole solutions 

The region of Kerr spacetime between the inner and outer horizons may be written as a 
Gowdy spacetime over S 2 x S 1 , by choosing as time coordinate a rescaling of the usual radial 
coordinate [ES|. In our parametrisation (^j, the velocity v tends to —1 along the axis and 



to +1 elsewhere. The discontinuous behaviour is an artefact of the parametrisation, since 
in the parametrisation (||]), DP — > —1 everywhere. This provides an example of Gowdy 
spacetimes on S 2 x S 1 without curvature singularities. The extremal case corresponding 
to the interior of the Schwarzschild black hole is a polarised Gowdy solution with v — > — 1 
at the horizon and v —>■ 3 at the curvature singularity. 

The same idea can be applied to the 'distorted black hole' spacetimes studied by Geroch 
and Hartle. These are constructed by analytic continuation of certain Weyl solutions, es- 
sentially corresponding to perturbations on an exterior Schwarzschild background, through 
the horizon. The solutions are static, but they can be different from Schwarzschild space- 
time by violating asymptotic flatness or by allowing matter in the exterior region. The 
interior is a polarised Gowdy spacetime with S 2 x S 1 topology. If we let the Schwarzschild 
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solution be given in Gowdy coordinates with P = P$, the Geroch-Hartle solution is given 
by P = Ps + U, where U is regular on the horizon. It is possible to show that U must be 
regular at the other singularity as well. Thus the asymptotic velocity is the same as for 
the interior Schwarzschild solution. 

Note that for polarised Gowdy spacetimes, analyticity is a necessary as well as sufficient 



condition for the existence of an extension with a compact Cauchy horizon [0. The 
necessity follows from the fact that on the extension, the evolution equation becomes the 
Laplace equation in suitably chosen coordinates. 

2.4.3 Numerical solutions 

Garfinkle has done numerical simulations of Gowdy spacetimes on S 2 x S 1 , with similar 
results as in the T 3 case [^]. All solutions have velocity —1 at the axis, without spiky fea- 



tures, which probably can be explained by the particular choice of initial data. We should 
emphasise however that Garfinkle's calculations are carried out in the parametrisation (Q), 
and rewriting his results in our parametrisation (|9]) gives solutions with 'false spikes' at 
the axis. Since we cannot generate solutions with spikes by the method in this paper, our 
results cannot be compared directly with those in [EyJ. 



3 The analytic case 
3.1 Fuchsian systems 

Consider a system of partial differential equations on R n+1 , whose solutions are expected 
to have a singularity as t — > 0. The Fuchsian algorithm is based on the following idea: de- 
compose the unknown into a prescribed singular part, depending on a number of arbitrary 
functions, and a regular part w. If the system can be rewritten as a Fuchsian system of 
the form 

td t w + N(x)w = tf(t,x,w,w x ), (23) 

where w x denotes the collection of spatial derivatives of w, the Fuchsian method applies. 
In the analytic case we have the following theorem. 

Theorem 3.1 (Kichenassamy and Rendall H28|| ). Assume that N is an analytic ma- 
trix near x = such that there is a constant C with \\o~ N \\ < C for < a < 1, where a N 
is the matrix exponential of Nina. Also, suppose that f is a locally Lip schitz function of 
w and w x which preserves analyticity in x and continuity in t. Then the Fuchsian system 



(^3|) has a unique solution in a neighbourhood of x = and t = which is analytic in x 
and continuous in t, and tends to as t — > 0. 

The method of proof is a variation of the standard Cauchy problem. Note that Theo- 



rem 



3TT| also applies to the case 

td t w + N(x)w = t a f(t,x,w,w x ), (24) 



since changing the t variable to t a transforms the system into the form 
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3.2 Low velocity, analytic case 

Let k, (p, loq and ip be real analytic functions of 9 on a neighbourhood of 9 = in S 2 , with 
k G (0, 1). We introduce new unknowns u and v such that 

Z(t, 9) = k(9) lnr + <p(9) + r e u(r, 9), (25a) 
u(t, 9) = u (9) + r 2fc W + v(r, 0)), (25b) 



where e > is a small constant which we will fix later. The idea is to use Theorem 3.1 



to show that it, v, Du and Dv tend to as t — > 0, so that ( p5[) provides an asymptotic 
expansion of Z and cj. The term 'low velocity' is then justified by the fact that the velocity 
v tends to 1 — k e (0, 1) as t — > 0, as is easily verified from the definition (|T0|). For this 
reason we will refer to 1 — k as the 'asymptotic velocity'. However, remember that from 
the discussion in section p73] , the velocity should be —1 or 3 at the axis, so the solutions 
constructed in this section are not regular at the axis. 
In terms of X, the expansion ( |25b| ) of uj corresponds to 

X = X {6) + t 2 ^-^ (4,(9) + v(r,9)), (26) 
for some functions X Q , %p and v, which can be compared directly with the T 3 results [28 



Inserting the expansion ( p5] ) into the system ( Jl6| ) gives 

(1 - t 2 )D 2 u = - (1 - r 2 ) [2eDu + e 2 u] + r 2 ' e [k + (Inr)AJfe + A<p] + r 2 [(D + e)u + Au] 

- e- 2 *- 2Ttu {(l - r 2 )r 2k -%D + 2k)(iP +v)] 2 - r 2 ~ 2k -\V uo ) 2 

(27a) 

- 2r 2 - e Vu ■ + v) + 2(lnr)(V> + v)Vk] 

- r 2 + 2k -t [v(V» + v) + 2Qnr)(ip + v)Vk]' 

(1 - t 2 )D 2 v = - 2k{l - t 2 )Dv + 2(1 - T 2 )r e (D + e)u(D + 2As) (-0 + v) 
+ t 2 [2(lnr)VA; • V(^ + v) + 2(Inr)(V> + v)Ak 

+ (D + 2k) (ij; + v) + A(4j + v)] (27b) 

+ r 2 - 2k Au - 2T 2 ~ 2k Vu ■ [(In r) Vk + V<p + r e Vit] 

- 2r 2 [W + r e Vu] ■ [2{\iit)(4j + v)Vk + + v)]. 

Next, we introduce the variables 

Uq := u, U\ := Du, U 2 := tu x , U 3 := rity, (28a) 
V :=v, V 1 :=Dv, V 2 := rv x , V 3 := rv y . (28b) 

We also write U and V as shorthands for the vectors (U 2 ,U 3 ) and (V 2 , V 3 ). Using fll9f) and 
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( pOD the system (p7[) may then be written as 



DU = U u (29a) 
(1 - t 2 )DU x = - (1 - r 2 ) [2eC/i + e 2 f/ ] + r 2 ~ e [k + (\nr)Ak + Atp] + r 2 (?7 a + eC/ ) 
+ r [(1 - x 2 )£/ 2:E - xy(C/ 22/ + £/ to ) + (1 - y 2 )U 3y - 2xU 2 - 2yU 3 ] 

- e-*r-*- u °{(l - rV k ^[V x + 2£# + V )] 2 - r^ k -\Vu,f 

- 2r 1 - £ Vcu ■ \V + rV^ + 2r(lnr)(^ + F )VA;] 

- r 2fc " e [F + rVf + 2r(lnr)(^ + V )Vk] 2 j, 

TO = r([/ to + C/ico), (29c) 

DU 3 = r(U 0y + U ly ), (29d) 

= 7i, (29e) 
(1 - t 2 )DV 1 = - 2k(l - r 2 )V x + 2(1 - T 2 )r e {U x + eU ) (Vx + 2k{^ + V )) 

+ t 2 [Vx + 2k(tP + V Q ) + 2(lnr)(^ + V )Ak + 2(\nr)Vk ■ + Aip] 
+ 2r(\nr)Vk-V + r 2 - 2k Au;o 

+ r [(i - x 2 )^. - xy(y 2 , + v 3x ) + (i - y 2 )y 3y - 2xv 2 - 2yv 3 ] u 

- 2T 2 - 2k Vuj ■ [(lnr)Vife + Vy? + r e W ] 
-2r[rV^ + r e Z7] • [2(lnr)(V> + V )Vk + V(V> + V )], 

DV 2 = t{V 0x + V 1x ), (29g) 
W 3 = r(yo» + n tf )- (29h) 

If we choose e such that < e < min{2A;, 2 — 2k}, the system ( |2"9"D is of the form 
(p4j) for some a (after dividing (|29b|) and ( ]291| ) by 1 — r 2 ), and a rescaling of r gives a 
Fuchsian system of the form (P3|) . By the criterion in IB], or by direct computation of a N , 
the conditions of Theorem |3J] are fulfilled. (Note that the essential features of the system, 
in particular the matrix N and the powers of r, are the same as in the T 3 case |2H[| .) We 
come to the following conclusion. 

Theorem 3.2. Suppose that k, if, uq and ip are real analytic functions of 8 in a neigh- 
bourhood of 9 = in S 2 such that k e (0, 1) for each 9. Let e be a positive constant less 
than min{2A;, 2 — 2k}. There exists a unique solution of Einstein's equations ([TT] ) of the 
form (f25|) in a neighbourhood of 9 = and t = such that u, v, Du and Dv all tend to 
as t — > 0. 



3.3 Negative velocity, analytic case 

In the case when we only assume k to be positive, allowing values of k greater than 1, 
we may use the same expansion fl2~5|) as in the low velocity case. However, examining the 
system (|2D|) we see that we must assume that ujq is constant in order to avoid negative 
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powers of r. This is again similar to the T 3 case, and is to be expected since numerical 
simulations indicate that spiky features may appear as t — > |28|, p~5]| . Again, using the 
definition (|10|) of the velocity shows that the asymptotic velocity is 1 — k, which is negative 
if k > 1. The expansion (|25b|) of lu implies 

X = X (9) + T e v(r,9). (30) 



for some Xq and v, which is in agreement with [^| . 



By the same arguments as in the low velocity case we have a similar existence theorem 
for a smaller set of data, including the regular case with velocity —1 at the axis. 

Theorem 3.3. Suppose that k, <p andip are real analytic functions of 8 in a neighbourhood 
of 9 = in S 2 such that k > for each 9, and suppose that uj is a constant. Let e be a 
positive constant less than 2k. There exists a unique solution of Einstein's equations ( |TTD 
of the form in a neighbourhood ofd = and t = such that u, v, Du and Dv all tend 
to as t — > 0. 

3.4 High velocity, analytic case 

If k is negative, we replace the expansion fl2"5|) with 



Z(t, 9) = k{9) In t + (p(6) + r e u(r, 9), (31a) 
uj(t,9) =oo (9) +r e v(r,9). (31b) 

Note that we do not get the full number of free functions in this case either. Again, this is 
because of the possibility of spiky features in the asymptotic behaviour. The calculations 
are similar to the k > case and will not be repeated here. It turns out that the system is 
valid for k < 1/2 (corresponding to an asymptotic velocity greater than 1/2), by choosing 
max{0, 2k} < e < min{2, 2 — 2k}. The corresponding expansion of X is 

X = X + T 2 ^ k ^^(9)+d(r,9)), 

for some functions ip and v and an integration constant X (which must vanish for the 
solution to be smooth at the axis). These solutions include the regular case with velocity 
3 at the axis. 

Theorem 3.4. Suppose that k, ip andu are real analytic functions of 9 in a neighbourhood 
of 9 = in S 2 , such that k < 1/2 for each 9. Let e be a positive constant such that 
max{0,2/;;} < e < min{2, 2 — 2k}. There exists a unique solution of Einstein's equations 
(|TT|) of the form (|3l|) in a neighbourhood of 9 = and t = such that u, v, Du and Dv 
all tend to as t — > 0. 
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4 The smooth case 



In the previous section we established existence of real analytic solutions with the desired 
asymptotic behaviour in a number of cases. Here we will consider the corresponding smooth 
solutions, using a generalisation of the approximation scheme in |38| . 

The basic idea is to approximate smooth asymptotic data (k, ip, u , if)) by a sequence 
of analytic data (k m ,<p m ,uJo m ,if> m ), and show convergence of the corresponding analytic 
solutions w m to a smooth solution w. Since the argument does not depend on the details 
of the Gowdy models, we will study a more general symmetric hyperbolic system with a 
singular term and sufficiently well-behaved coefficients. Also, since we are only interested 
in a neighbourhood of one of the axes in S 2 x S 1 or S 3 , we will carry out the argument on 
a subset of [0, oo) x R n . 

Definition 4.1. We say that the system of differential equations 

tA°(t,x)d t w + N(x)w + tA j (t,x,w)d j w = tf(t,x,w) (32) 

is regular symmetric hyperbolic if A is uniformly positive definite and symmetric, the A^ 
are symmetric, and all coefficients are assumed to be regular in a sense to be specified 
below. 



Note that our terminology differs from the usual one here, since the full system ( j32|) is not 
regular at t = in the usual sense. In [[38|] , the case when A is the identity was considered. 
We will therefore only give an outline of how to generalise the argument to the case of more 
general A . 

In our case, the system (^) can be written as the Fuchsian system 

td t w + N(x)w = tf(t,x,w,w x ). (33) 

Indeed, if A were independent of t, we could put N := (A°)~ 1 N and / := (A°)~ 1 (f — 
A^djiv). Now the t-dependence of A in our case is such that we may write (A°) _1 (t, x) = 
B (x) +t a Bi(t,x) for some regular matrices B and B\ and a constant a for small enough 
t. The t dependence of (t4 ) _1 A^ may then be included in / (after rescaling t if a < 1). 



The techniques of |38[ can then be applied to (£33|) to obtain formal solutions (see below). 

Note also that most of the arguments below hold in the more general case when A 
depends on w as well. However, in that case it is often not possible to rewrite the symmetric 
hyperbolic system Fuchsian system of the form 



4.1 Regularity and formal solutions 

Here we recall some results from |38| on the existence of formal solutions. We first need 
appropriate notions of regularity and of formal solutions. 

Definition 4.2. A function f(t,x) from an open subset Q C [0, oo) x R n to R m is called 
regular if it is C°° for all t > and if its partial derivatives of any order with respect to 
x G R n extend continuously to t — (within Q). 



14 



Definition 4.3. A finite sequence (wi, W2, ■ ■ ■ , w p ) of functions defined on an open subset 
Q C [0, oo) x R n is called a formal solution of order p of the differential equation (|33|) on 

n if 

1. each Wi is regular and 

2. td t uii + N(x)wi — tf(t, x, Wi, Wi X ) = 0{t l ) for all i as t — > in f2. 



The existence of formal solutions to a system of the form ( j33|) is provided by the 
following result. 



Lemma 4.1 (Rendall |]38[j). /// is regular and N is smooth and satisfies \\c N \\ < C for 
some constant C and all a in a neighbourhood ofO, then ( |33"D has a formal solution of any 
given order which vanishes at t = 0. 

By the discussion in the previous section, a formal solution of (Bl) may also be regarded 



as a formal solution of ( p2|) with a definition similar to Definition [4.3| . We have thus 



established the existence of formal solutions to ( p2|) of any order. 



In proving the existence of smooth solutions to ([32] ) it is important that the matrix N 
is positive definite. This is not the case in general, so we need to modify the system to 
fulfil this requirement. The idea is to use a formal solution {wi, . . . , Wj} and to study the 
system satisfied by Z{ := t l ~ l {w — w^, where w is the sought solution of the original system 
(|32|). The procedure is similar to that in |ft8| , so the details will be omitted. We end up 



with a system 

tA°(t, x)d tZi + (N{x) + l)A°{t, x)) Zi + tA j (t, x, Wi + f'hJdjZi = tfi(t, x, Zi ), (34) 

where the regular function /j is constructed from and / and depend on Wi and Wi X as 
well. The point is that since A is uniformly positive definite, the coefficient iV+ (i — 1)A° 
is positive definite for large enough i. We will refer to fl3~4]) as the positive definite system. 

4.2 The existence theorem 

Given smooth asymptotic data S := {k, <p, u , ip} on U C R n , we may construct a sequence 
of analytic data S m := {k m , ip m , uj 0m , ip m } on U which converges to S in C°°(U), uniformly 
on compact subsets. If the formal solutions are constructed as in the proof of Lemma |4.1| 
(see |^8|), the analytic formal solutions w m i of order i corresponding to the analytic data 
S m converge to a formal solution Wi of order i corresponding to the smooth data S as 
m — > oo, and the convergence is uniform on compact subsets. This also holds for spatial 
derivatives of any order. Hence spatial derivatives of any order of the coefficients of the 



positive definite system (|34|) converge on compact subsets as m — > oo. It follows that on 
any compact subset there is an i such that the coefficient involving N(x) is positive definite 
for all m in our case, so we fix such a value of i and omit the index i from now on. 

The global existence theorem for Gowdy spacetimes over S 2 x S 1 and S 3 |]TU| implies 
that there are smooth solutions of (153) on a common time interval for all m. Thus our 



problem can be solved by proving the following theorem. 
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Theorem 4.1. Let z m (t, x) be a sequence of regular solutions on [0, ti) X U C [0, oo) x R n , 
with z m (0,x) = 0, to a sequence of regular symmetric hyperbolic equations 



^m(^' x )dt z m + N m (x)z m + tA^t, X, Z m )djZ m — tf m {t, x, z m ). 



(35) 



Suppose that N m is positive definite for each m and that the coefficients converge uniformly 
to Aq, N , and fo on compact subsets as m — > 0, with the same properties as A m , N m , 
A J m and f m , and that the corresponding spatial derivatives converge uniformly as well. Then 
z m converges to a regular solution z of the corresponding system with coefficients Aq, N , 
Aq and f on [0, t ) x U for some t , and z (0, x) — 0. 

Proof. The idea is to use energy estimates to show that {z m } is a Cauchy sequence. Since 
the proof is very similar to the case when A is the identity []38|1 , we only give an outline 
of the important steps here, with emphasis on the differences. 

First we consider the system satisfied by spatial derivatives of z. The collection z of 
spatial derivatives up to order s satisfies a system similar to (|35j), 



^■^■rrSti x )^t z m + N m (x)z m + tA^ )1 (t, X, Z m )dj"L m — t f m (t, X, Z m ) , (36) 

obtained by differentiating fl3"5p and substituting the equations for lower order spatial 
derivatives. The problem is that N m is not necessarily positive definite, due to the presence 
of off-diagonal blocks depending on (A^J -1 and spatial derivatives of A° m and N m . But 
this can be dealt with by multiplying the spatial derivatives D a z m by for a sufficiently 
small constant K as in |$8| (here a is a multi-index). Let z m := {K^ D a z m ; \a\ < s} be 
the collection of weighted spatial derivatives up to order s and assume that A m , N m , 
and f m are the corresponding weighted coefficient matrices. For example, in one spatial 



dimension the system for z r 



d T z m ) has coefficients 



A° m 
A® 



N m 
KD m N m N n 



Al 
KD m Al Al 



f 



fm 

KD m f n 



(37) 



where D m is the operator d x — (c^A^)^^) -1 . The convergence of the coefficients for any 
fixed s follows from the convergence of the coefficients of the original equation. 
We will need appropriately weighted Sobolev norms 



\u 



IE 

\a\<s 



(38) 



Note that these are equivalent to the usual Sobolev norms. The norms depend a priori on 
m since they include a factor A^, but since A^ converges uniformly on compact subsets 
to A which is uniformly positive definite, the equivalence can be taken to be independent 
of m on compact subsets for sufficiently large m. 

Second, a domain of dependence result may be obtained by standard techniques. The 
argument is the same as that in |38| and will be omitted. Thus we need only consider the 
problem on a compact subset. 
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Next, we will show that the sequence ||z m ||#s is bounded. Differentiating and using the 
definitions of z m and we get 

d t (\\z m \\ 2 Ha ) = 2(A° m d t z m ,z m ) L 2 + ^ K 2 ^(D a z m , (d t A° m )D a z m ) L 2, (39) 

\a\<s 

and using fl3~6f) to substitute for A^<9 t z m gives 

dt{\\z m \\ 2 H s) = -2r 1 (N m z m ,z m ) i2 + R m , (40) 

where R m contains the same terms as in the regular case with N m = 0. The first term 
on the right is negative since N m is positive definite by construction, and R m may be 
estimated as in the regular case, giving 

H- Call Zm || H s ■ 

(41) 

Here d and C 2 are polynomials in || (vl )— 1 1| ^oo , \\D a A°\\ L oo, \\D a A^\\ L oo and \\D a f\\ L °o for 
|a| < s, and in ||z TO i|#/s for any given k > n/2 + 1 (by the Sobolev embedding theorem). 
Since all of the coefficients converge and we may choose k = s if s > n/2 + 1, dt(||2: m ||fp) 
is bounded by a polynomial in ||z m ||/p whose coefficients are independent of m. Applying 
Gronwalls lemma then gives that ||£ TO ||fl" s is bounded. 

We can in fact obtain a stronger estimate. Since ||z m ||#s is bounded, it follows that 
the coefficients C\ and C 2 are bounded, so they can be chosen to be constants. A second 
application of Gronwalls lemma then gives that £ _1 ||z m ||i2" s is bounded if s > n/2 + 1. 

Finally, we show that {z m } is a Cauchy sequence in the H s norm. Let the difference 
between consecutive elements of the sequence be given by v m := z m — z m -\ and put 
v m := z m - z m _i. From @ we get 

|^(lkm||| s ) = - t~ l (N m v m , V m ) L 2 - t _1 ((N m - N m _i)z m _i, V m ) L 2 

- (( A m ~ Kn-X^m-U Vm) V* ~ (A-L( Z m) d jV m , V m ) L 2 (42) 
+ (v m ,g m ) L 2, 

where 

g m := t(i m {z m ) - f m _i(z m _i)) - t(A J m (z m ) - A? m _ 1 {z m ^ 1 ))d j 2 lm _ 1 . (43) 

The first term on the right hand side of (|42"D is negative and can be discarded, and using 
that t -1 ||,2 m _i||H s is bounded, the second term is less than C||N m+ i — N m ||z°°||'Urn||ff s - To 
estimate the third term we need a bound on ||£^z m _i ||# s , but this can be obtained by 
applying the bound on £ -1 ||2m-i||.ff«+i to (|55|). The fourth term may be estimated by doing 
a partial integration as in the regular case. The last term is estimated by inserting some 
terms and applying the mean value theorem and the bound on ||z m _i||^ s +i. We end up 
with an estimate of the form 

^(IkmlliP-) < Ci|| -h Call v m \\H", (44) 
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where C\ is a constant independent of m and C 2 is a polynomial in ||A m — A m _^||^oo, 
||N m — N m _i ||x,«>, || A J m — A^^j^H^oo and ||f m — f m -i which tends to as m — > 00. 
Applying Gronwalls Lemma again, we conclude that ||u TO ||/rs — > as m — > 00. 

We have shown that the solutions z m of the positive definite system (|35|) converges in 
the .fP norm to a smooth solution for each s. A potential problem is that the time interval 
of convergence may depend on s. But by standard techniques for symmetric hyperbolic 
equations the solutions may be extended to a common time interval, so we can conclude 
that z m converges as smooth functions as well. □ 



4.3 Low velocity, smooth case 

From the discussions in the previous section it follows that we only need to write the Gowdy 
equations ( p9|) in symmetric hyperbolic form, since the existence of a smooth solution is 
then guaranteed by Theorem We rewrite Q2"5|) as 

DU Q = U u (45a) 
(1 - r 2 )DU x = - (1 - r 2 ) [2el7i + e 2 U ] + r 2 " e [k + (lnr)AA; + Aip] + r 2 ^ + eU ) 
+ r [(1 - x 2 )U 2x - xy(U 2y + U 3x ) + (1 - y 2 )U 3y - 2xU 2 - 2yU 3 ] 

- 2 *- 2 ^{(l - r 2 )^-^ + 2A# + Vo)} 2 - r 2 - 2k -^Vu ) 2 

- Tt^Vuq ■ \V + rVf + 2r (In r) (V> + V )Vk] 

- r 2k - e [F + rV^ + 2r(ln r) (ip + V )Vk] 2 

(1 - x 2 )U 2 - xyU 3 + r(l - x 2 )U lx - rxyU ly , (45c) 
(1 - y 2 )U 3 - xyU 2 + r(l - y 2 )^, - rxyU lx , (45d) 

(45e) 



— e 



1 - a: 2 )/)^ 


- xyDU 3 


1 - y 2 )DU 3 


- xyDU 2 


DV 


= v u 


1 - T 2 ) 


= - 2k{l 



+ t 2 [Vi + 2k(^P + V ) + 2(lnr)(^ + Vb)Afc + 2(lnr)Vfc- + A 



+ 2r(lnr)Vfc ■ 1/ + r 2 ~ 2k Auj, 







(45f) 



+ r [(1 - x 2 )V 2x - xy(V 2 , + V 3 x) + (1 - y 2 )V 3y - 2xV 2 - 2yV 3 ] 

- 2r 2 - 2fc Vcu • [(In r)Vk + Vip] - 2r l+e - 2k Vu Q ■ U 

- 2r e fr 1- ' + U] ■ [2r(ln r) (V> + Vo)V/c + rV^ + 7] , 

(1 - x 2 )DV 2 - xyDV 3 = (1 - x 2 )V 2 - xyV 3 + r(l - rr 2 )^ - rxy^, (45g) 
(1 - ?/ 2 ) W 3 - a^W 2 = (1 - y 2 )V 3 - xyV 2 + r(l - y 2 )^, - rxyV lx . (45h) 

Because of the factor f 1+e ~ 2k , we must restrict k to the case < k < 3/4. After 
an appropriate rescaling of the time variable, the system is clearly a regular symmetric 
hyperbolic system of the form fl32|). Forming the corresponding positive definite system 
(where i = 3 suffices to make the coefficient of the singular term positive definite) and 
applying Theorem |4.1| then gives the existence of the desired solutions. 
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Note that there is a slight difficulty in showing that ( [45| ) is equivalent to (p9|). For 
example, from ( [45c|) , ([45d|) and ( |45a|) it follows that D(U 2 — Td x U ) — U 2 — Td x U , and 
since U2 and C/o vanish at r = we must have U% — tO x Uq = c(x)t for some function c(x). 
But there seems to be no direct way of showing that c(x) = follows from fl45|) . This 
must be the case however since we obtain the solution as the limit of a sequence of analytic 
solutions with this property 

We summarise this section in the following theorem. 

Theorem 4.2. If k, ip, ujq and ip are smooth functions of 9 and < k < 3/4 for all 9, 
there exists a solution of Einstein's equations (|Tl"D in a neighbourhood U of 9 = of the 
form fl25|) where 2k — 1 < e < min{2A;, 2 — 2k} and u and v are regular and tend to as 
t — > 0. Given the form of the expansion and a choice of e, the solution is unique. 



4.4 Negative velocity, smooth case 

Since exponents of r involving —2k only appear in terms containing Vuj , the argument 
in the previous section applies immediately to the negative velocity case with k > and 
constant c^o- 

Theorem 4.3. If k, ip and ip are smooth functions of 9 such that k > for all 9 and lvo 
is a constant, there exists a solution of Einstein's equations in a neighbourhood U of 
9 = of the form (|25D where e < 2k and u and v are regular and tend to as t — + 0. Given 
the form of the expansion and a choice of e, the solution is unique. 

4.5 High velocity, smooth case 

The calculations for the high velocity case are similar to the low velocity case and will be 
omitted. We have the following existence result. 

Theorem 4.4. If k, ip anduo® are smooth functions of 9 and k < 1/2 for all 9, there exists 
a solution of Einstein's equations (|ll|) in a neighbourhood U of 9 = of the form ([31]) 
where max{0, 2k} < e < min{2, 2 — 2k} and u and v are regular and tend to as t — > 0. 
Given the form of the expansion and a choice of e, the solution is unique. 

4.6 Intermediate velocity, smooth case 

It remains to treat the case when 3/4 < k < 1. The idea is to include one more term in 



the expansion (|25a| ) of Z, 

Z(t, 9) = k(9) lnr + (p(9) + a(9)r 2 - 2k + r 2 - 2k+t U (r, 9), (46) 

where a is to be chosen such as to eliminate the leading order terms in the equations. We 
keep the original expansion for u as in ( |25b ). 
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The resulting system is similar to ([45]), but with DU\ and DU 2 given by 



r 2 )DU 1 



- (1 - r 2 ) [(4 - 4k + 2e)U x + {2-2k + e) 2 U Q ] - 4r(ln r) Vk ■ U 
+ r 2k ~ e [k + (In r)Ak + Aip] + T 2 - £ {4k 2 - 10k + 6)a 
+ r 2 - £ [Aa - 2(lnr)(aAk + 2Va • Vk) + 4(lnr) 2 a(VA;) 2 ] 
+ t 2 [Ui + (2-2k + e)U - 2(\nr)U Ak + 4(lnr) 2 f/ (VA;) 2 ] 
+ r[(l - x 2 )U 2x - xy(U 2y + U 3x ) + (1 - y 2 )U 3y - 2xU 2 - 2yU 3 ] 



exp(— 2y? 
'(1 

2r 



2ar 



2-2k 



2r 2 - 2k+e U ) 



(47a) 



ik-2-e [ 

T IT 



+ T~ 



[Vi + 2fc(</> + Vb 
2fc - 1 - e Va;o ■ [F + tV^ + 2r(lnr)(V> + V )Vk] 
- r 4fc - 2 - e [F + r v^ + 2r(ln r) (V> + V ) Vk] 
[exp(-2^ - 2«T 2 - 2fe - 2r 2 - 2k+e U )(Vuj ) 2 - 
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■ [2r(\nr)(ip + V )Vk + rVip + V], 

We choose a := (2 — 2&;)~ 2 (Va;o) 2 exp(— 2yj) to cancel the T _e factor in the last term 
on the right of ( |47a| ). That term will then be of order r 2 ~ 2k ~ e , S o the order of the (Vuq) 2 
term is unchanged from the previous case ( [45b ). 

The problematic exponents of r are now 2 — 2k — e, 3 — 4k + e and 2k — 1 — e. We have 
to choose the number e > such that 4k — 3 < e < 2k — 1, but we also need to ensure that 
2 — 2k — e > 0, which is compatible with 4k — 3 < e if and only if 1/2 < k < 5/6. In that 
case, we may apply the arguments of section |44] and [O] to show that a regular solution 
exists. 



Contrary to what was claimed in ||38|| , it is not possible to cover the whole range 
1/2 < k < 1 since we can only choose a such that the last term in ( |47a| ) vanishes to first 
order in r. In equation (27) of |38| the corresponding higher order terms have been left 
out. 

We can however cover small intervals of k closer to 1 by repeating the method above. 
Replacing the expansion ([25 a|) by ([46]) is equivalent to performing the transformation 
u 1— > ar 2 ~ 2k ~ e + r 2 ~ 2k u. Repeating this transformation % times, where each a is chosen 



20 



such that the leading order coefficient of (Va>o) 2 is cancelled at each stage, gives a system 
with positive powers of r if and only if 



1 - — < k < 1 — (48) 

2i 2i + 4 ( J 

and 

1 - 2(i + 1)(1 - k) < e < min{2 - 2k, 1 - 2i(l - k)}. (49) 

The interval (1/2, 1) is covered by the infinite sequence of intervals (1 — 1/2?, 1 — l/(2z+4)). 
The expansion of Z for a given % is 

i 

Z(t, 9) = k(9) lnr + ip{9) + ^ a,(0)<r (2 - 2 ^' + T^- 2k ^ i+e u(r, 9). (50) 

j=i 

Note that as tends to 1, we have to include an increasing number of terms in the 
expansion ( |50|) of Z. This might seem to be contradictory since we keep the original 
expansion of u, and the higher order terms of uj should affect Z at some finite order. But 
for any given k, the number of terms in ( |50| ) is finite since % is bounded above by (|4*8|). 
Also, the exponent of r in the terms containing <x,- is always between and 1, so higher 
order contributions are still encoded in u. 

The above discussion motivates the following theorem. 

Theorem 4.5. If k, (p, u and ip are smooth functions of 9 and k satisfies ([Rj) for all 
9 and some natural number i, there exists a solution of Einstein's equations (|TT|) in a 
neighbourhood U of 9 = of the form ( |50D and (|25b|) where e satisfies (^9|) and u and v 



are regular and tend to as t — > 0. Given the form of the expansion and a choice of e, the 
solution is unique. 



5 Discussion 

We have constructed families of solutions to the evolution equations of Gowdy spacetimes 
with S 2 x S 1 or S 3 spatial topology. When the asymptotic velocity is between and 1 
(0 < k < 1) we obtain solutions depending on four free functions, which is the same 
number as in the general solution, while outside that range the solutions include only 
three free functions. The solutions are asymptotically velocity dominated by construction. 
Unfortunately, for regular asymptotically velocity dominated solutions of the full Einstein 
equations the asymptotic velocity must be —1 or 3 (k = ±2) at the axes. This can also be 
checked directly by inserting the asymptotic expansions of the solutions into the constraints 
(|17|). Since both analytical and numerical arguments indicate that the velocity is between 
and 1 for T 3 Gowdy, and these results may be applied on sets not containing the axes, 
it seems that there will be spikes at the axes in general situations. 
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There are two possible interpretations. Firstly, there is the possibility of false spikes, 
which result from an inappropriate choice of parametrisation of the metric. This is also 



what makes our solutions different from the numerical solutions [15]. The parametrisation 
problem may be traced to different choices of parametrisations of the hyperbolic plane. It 
might be possible to choose this parametrisation differently such as to avoid coordinate 
singularities. This is of equal importance in the study of the dynamics of T 3 Gowdy, in 
particular when trying to verify the numerical observation that the velocity is eventually 
driven below 1 in a mathematically rigorous way. 



Secondly, there is the possibility of true spikes at the axes. In [39], T 3 Gowdy solutions 
with spikes are constructed by applying suitable transformations to a given solution with 
velocity between and 1. At first sight it seems possible to do something similar for the 
other topologies. The transformations used in p9| are inversion, which interchanges the 
Killing vectors, and a Gowdy-to-Ernst transformation. A combination of these transfor- 
mations preserve the evolution equations in the S 2 x S 1 and S 3 C9..SCS clS well. In fact, the 
equations for Z and u are the same as the equations for — P and Q (in the notation of 
section |2.1| ). But P is singular at the axes, so this method of constructing solutions with 
spikes does not work in our case. 

Finally, some remarks on the asymptotic behaviour of the curvature is in order. Using 
the asymptotic expansions for our solutions, (|25|) or (pi]), it is straightforward to calculate 
that a necessary condition for the Kretschmann scalar RijkiR 1 ^ 1 to be bounded is that the 
asymptotic velocity is ±1 (k = or k = 2), which is in agreement with the result for the 
polarised models described in section [2.4. 1| . All of the black hole solutions mentioned in 



section |2.4.2j are of course extendible through the horizon, and the asymptotic velocity is 
indeed ±1 there. While such non-generic solutions can be ignored from the cosmological 
point of view, they might be interesting by analogue to black hole interiors. One interesting 
question that remains to be answered is under which circumstances general Gowdy models 
can be extended through a compact Cauchy horizon. In particular, is analyticity necessary 

H? 



as in the polarised case [JT2 



To conclude, it seems that for more progress on Gowdy S 2 x S 1 and S 3 spacetimes to 
be made, a better understanding and handling of the spikes are needed. If this can be done 
for T 3 Gowdy, it should be possible to adapt the techniques to the other topologies using 
some of the arguments of this paper. 
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